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I. INTRODUCTION 



There are two inspiring notions in quantum mechanics called supersymmetry and superin- 
tegrability. Being formally independent, both of them are guide signs in searching for exactly 
solvable problems. Moreover, some of quantum mechanical systems, like the Hydrogen atom or 
isotropic harmonic oscillator, are both superintegrable and supersymmetric. Let us note that 
just such systems as a rule are very interesting and important. 

A quantum mechanical system with n degrees of freedom is called superintegrable if it admits 
more than n — 1 integrals of motion. The system is treated as supersymmetric in two cases: 
when some of its integrals of motion form a superalgebra, and when its Hamiltonian has a 
specific symmetry with respect to the Darboux transform, called shape invariance 

The search for superintegrable systems was started with paper where all the linear and 
quadratic integrals of motion for the 2d Schrodinger equation with an arbitrary potential had 
been presented. The systematic search for systems whose integrals of motion are first and 
second order polynomials in momenta was performed in B, B and Q. 

In last decades a number of new important results in this field has been obtained. In 
particular, integrable systems with second and third order integrals of motion in 2d and 3d 
Euclidean spaces had been intensively studied , the systems both electric and magnetic 
external fields and integrable systems with spin were discussed in ly] and HI], 12|, 13|. 

However, there exist interesting integrable systems which were not studied systematically 
till now. They are neutral particles with non-trivial spin and dipole moment (e.g., neutrons) 
interacting with an external electromagnetic field. These particles are described by Schrodinger- 
Pauli equations which include the Stern-Gerlach or electric dipole terms. A perfect example 
of such system is the Pronko-Stroganov model [l^ describing a neutron coupled to the field 
of the constant straight line current. This 2d system admits four integrals of motion (includ- 



ing HamiltonianJ 
the relativistic 



14| . three of which are algebraically independent. Notice that there exist 



15| and arbitrary spin ^] versions of this system which are integrable and 
supersymmetric too fl7 |. 

In addition, the Pronko-Stroganov system is shape invariant and can be easily integrated 
using tools of supersymmetric quantum mechanics 18|], 19|. This circumstance had motivated 
us to search for other supersymmetric Schrodinger-Pauli equations and classify matrix shape 
invariant potentials. In this way a number of new exactly solvable systems had been found j20| , 
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2l|. 

In the present paper the integrable 2d Schrodinger-Pauh equations for neutral particles 
are classified. Like in paper ll| we restrict ourselves to the first order integrals of motion 
and classify all external fields which give rise to integrable and superintegrable systems. As 
a result a a new class of integrable systems has been found. The majority of these systems 
is supersymmetric also, and one of them is shape invariant. Thus the germaneness between 
integrability and supersymmetry becomes apparent also in quantum mechanical models of 
neutral particles. 

Integrals of motion of integrable systems presents powerful tools for finding their exact 
solutions. In this paper we restrict ourselves to solving three of the obtained systems. One of 
them appears to be a special case of models with shape invariant effective potentials that have 
recently been classified in 20|| and 2l|. The other system describes the neutron interacting 
with a periodic external field, which have both discrete and band energy spectra. One more 
system is rotationally invariant and includes a superposition of two external fields. 



II. CLASSIFICATION PROBLEM 



Let us consider a special class of Schrodinger-Pauli equations describing neutral fermions 
with non-trivial dipole momentum interacting with an external field. The corresponding sta- 
tionary Schrodinger-Pauli equation looks as follows: 

if^(x) = Eip{x.) (1) 

where 

H=(^ + ^cr.B]. (2) 

Here cr is the matrix vector whose components are Pauli matrices, B = B(x) is a vector of 
magnetic field strength, and vector x represents independent variables. In addition, A denotes 
the constant of anomalous coupling which is usually represented a.s X = gfio where fio is the 
Bohr magneton and g is the Lande factor. 

We shall classify integrable systems ([T]), i.e., find all Hamiltonians (E]) admitting a sufficient 
number of integrals of motion which are linear combinations of momenta with matrix coeffi- 
cients. In this paper we restrict ourselves to planar systems depending on two variables xi and 
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X2- However, bearing in mind possible generalizations of the presented results to 3d systems, 
we will not exclude the third coordinate a priori. Till an appropriate moment our analysis will 
be valid for both 2d and 3d systems. 

To obtain more compact formulae let us rescall variables and reduce Hamiltonian ([2]) to the 
following form: 

H = -V^ + o- ■ B (3) 

where V is the gradient vector with componets Vq = jf-. To achieve this goal it is sufficient 

OX a 

to change in 1^ E ^ and B ^B. 

For 3d systems = ^ + ^ + ^ and B depends on three variables Xi, X2, and X3. In the 
case of 2d systems the third variable X3 and the corresponding derivatives should be deleted. 
This action can be formalized by imposing the constraint V3'?/'(x) = and looking for integrals 
of motion which commute with H and V3. 

Let us search for integrals of motion for hamiltonian ([3]) of the following generic form: 

Q = a^(i{A^^VJ + ^]^) (4) 

where summation is imposed over the repeated indices = 0, 1, 2, 3 and a = 1, 2, 3 or a = 1, 2 for 
3d or 2d systems respectively. A''" and are functions of x, {A'''', Va} = A'"'Va+VaA''", Va = 
0"^ are Pauli matrices: 



By definition, integrals of motion should commute with Hamiltonian, 

[H, Q] = HQ - QH = 0. (5) 

Substituting ([3]) and (jlj) into ([5]), using the relations 

{a", a''} = 25"^ K> a^] = 2ie"''V'=, a, 6, c = 1,2, 3 

and equating coefficients for linearly independent matrices and differential operators, we obtain 
the following system of determining equations for coefficients A'^" and fi^: 







Ar + K' = 0, (6) 



nl = 0, (7) 

^ab^a ^ ^Ofe^a _ ^abc^b ^ (g^, 

Ql = 26^"^A^^B'^. (9) 

Here the subindices denote derivatives w.r.t. the corresponding independent variables, i.e., 
= etc., and summation is imposed over the repeated indices. 

The system of equations (EJ-Q presents the necessary and sufficient conditions for commu- 
tativity of operators ([3]) and (j4]). Equations and ([7]) are easy integrated, and their solutions 
have the following form: 

were CM'^^l,C'^''and are arbitrary constants satisfying the condition C^^""^^ = — (7^'^'*] . Substi- 
tuting f lTOj) into dH]) and (Q we obtain the overdetermined system of ffist order partial differential 
equations for functions Q'^ and B'^. 



III. DETERMINING EQUATIONS AND EQUIVALENCE TRANSFORMATIONS 
FOR 2D SYSTEMS 

The determining equations (jH])^© and their partial solutions f llUp are valid for both 2d 
and 3d equations ([T]). Starting with this point we restrict ourselves to two dimension systems 
depending on variables xi and X2. The corresponding function A'^" in ( ITO|) and the related 
operator @ are reduced to the following forms 



^^,a ^ C^^ba^^ + 0^"^ (11) 

Q = a^ {C^L + C^'^Pa) + a'^fi" (12) 

where and C^" are arbitrary real constants, a = 1, 2, b = 1,2, /i = 0, 1, 2, 3, Pa = 
— iVa, L = X1P2 — X2P1, s""^ = —e^"" and = 1. Moreover, the number of arbitrary constants 
in (fTT]) and ( !T2|) can be reduced using the equivalence transformations which keep the form of 
Hamiltonian (jS]) up to multiplier p-: 

Xa^ Xa + Ca, (13) 
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Xa^RabXb, a, 6=1,2, (14) 

_^ RkuQu^ _^ ^kn^n^ k,n = l, 2, 3, (15) 

Xa^XXa, B'^^B" (16) 

where i?^^ and R^"' are planar and spatial rotation matrices correspondingly, Ca and A 7^ are 
real constants. 

Let C° = a 7^ then up to shifts (1131) we can set = = 0. Moreover, up to rotation 
transformations (fTSl) we can restrict ourselves to = b,C^ = C'^ = 0. Then, applying 
rotations (HM of variables Xi, X2 and, if necessary, rotations ( |T5l) with k,n = 1,2 we can reduce 
the remaining constants C"^ to C^^ = ci,C^^ = C2,C^^ = = where ci, C2 and 6 are 
arbitrary parameters. As a result we reduce functions A'^" in ( ITT]) to the following form: 

A^i = ax2, A°^ = -axi, A^^ = 6x2 + c?i, 

(17) 

A32 = + d2, A" = Ci, A22 = C2 

while all the other components of tensors A'^'^ are zeros. In ( [T71) we denote constants C^^ and 
as di and ^2 correspondingly. 
If one of parameters a or 6 (or both of them) are zero, we can use (fT3l) again and reduce the 
set of functions f|T7|) to the following non-trivial components: 

A^^ = bx2, A32 = -bxi, A" = ci, A22 = C2, A^^ = C3, A°2 = C4, by^O (18) 
A°^ = ax2, A°2 = -axi, A^^ = ci, A^^ = C2, a 7^ 0, (19) 
Aii = ci, A22 = C2, A0i = C3, A02 = c4. (20) 

The corresponding symmetry operators ( |T2l) look as follows: 

Q = aL + asbL + criCiPi + (T2C2P2 + o"afia, 7^ 0, (21) 
Q = asbL + (C3 + aiCi)Pi + (C4 + o-2C2)P2 + crA, b ^ 0, (22) 
Q = aL + (JiCiPi + (T2C2P2 + o-^ria, a 7^ 0, (23) 

Q={C3 + CJiCi)Pi + (C4 + a2C2)P2 + CTa^a (24) 

where L = X1P2 ~X2pi, Xi = Xi — ^, X2 = ^2 + ^, Pa = — i^f-- In addition, without loss 
of generality we can set a = 1 in fl2Tl) and fl23|) and 6 = 1 in fl22|) . 

Notice that for ci = C2 = C5 = operators ( l23l) and (l24l) are reduced to the following form: 



g = - ^as, (25) 



Q = C3P1 + C4P2 - 2^3, C3C4 = (26) 

where k, n, C3 and C4 are real parameters. Indeed, in accordance with (Q, in this case functions 
are reduced to constants, and hermitian matrix (Jafi" is diagonahzable. Just operators fl25|) 
and (l26ll represent Lie symmetries which could be admitted by equation ([1]). 

Equations (l2T])-( !26l) give representatives of the family of operators ( fT2l) defined up to equiv- 
alence transformations (|T3l) -( fT6l) . In order these operators to be integrals of motion for Hamil- 
tonian functions B^,B'^,B^ and have to satisfy equations ([S]) and which, in 

view of f|T7j) - fl20|) . take the following form: 

Ql = -2ciB^, Ql = 2c2B\ 

nl = -2(6x2 + di)B\ nl = 2{d2 - bxi)B\ 

nl = 2ciB^ + 2{bx2 + di)B\ 

nl = -2c2B^ + 2(6x1 - d2)B\ 

^ ^ (27) 

6 {xiBl - X2BI) = ciB\ + C2BI, 

a{xiBl - X2BI) + n^B^ - n^B^ = c^B\ + dBl, 

a{xiBl - X2BI) + 9}B^ - V?B^ = c^Bl + c^Bl, 

a{xiBl - X2BI) + n^B^ - n^B^ = c^iBf + c^Bl 

where parameters di, C3, C4, a and 6 satisfy the conditions (2C3 — CIC4. = 0; di = d2 = 
if a6 = and + 6^ > 0. 

Thus the problem of classification of superintegrable planar Schrodinger-Pauli equations 
admitting first order integrals of motion is reduced to finding the general solution of determin- 
ing equations ( |2711 . These equations should be solved consequently for all sets of parameters 
a, 6, di, (^2, Ci, C2, C3, C4 present in operators (pT]) - (!26|) . Thus for operator of type (pT!) we set in 
fl27j) C3 = C4 = 0, a 7^ 0, 6 7^ while the remaining parameters are arbitrary, etc. 

IV. LIE SYMMETRIES 

Let us start with an important class of symmetry operators which are generators of con- 
tinuous groups. They are the special cases of operators f|T2|) corresponding to A^" = for all 
/i 7^ 0. Up to equivalence all such operators are given by equations (1251) and (l26l) . To classify 
the external fields B = {Bi, B2, B^) such that Hamiltonian ([3]) commutes with operators f l25l) 
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and (I2SD it is sufficient to solve equations f E7|) with Qi = Q2 = b = di = d2 = 0,^3 = ^ or 
^3 = |, while parameters C3,C4,n,k and a should satisfy the conditions a{cl + C4) =0 and 
nk = 0. As a result we obtain the list of the external fields and the corresponding symmetries 
presented in Table 1. A constant external field is not presented here since the corresponding 
Hamiltonian can be reduced to a direct sum of two Hamiltonians with trivial interaction terms. 



Table 1. External fields with Lie symmetries 



No 


External field 


Symmetry operators 


Comments 


-1 

1. 


fii = cos{ke)fi{r) + sin(fc^)/2(r), 
= cos(/c6')/2(r) - sm(/c6')/i(r), 
= fsir) 


VI = -t^ + 2*^3 


V • B 7^ 


2. 


= ficos{kO)r B"^ = fism[kO)r ^, 
5 = fsir) 






3. 


B^ = ficos{e)Mr), 52 = /isin(^)/i(r), 
= /2(r) 


Qi = L + icrg 




4. 


i?^ = cos{6xi)fi{x2) + sm{5xi)f2{x2), 
52 = cos(5xi)/2(x2) - sin(5xi)/i(x2), 

5' = f3{x2) 


Q2 = Pl- |(T3 


V-B ^0 


5. 


B^ = yuexp(— X2) cosxi, 
-B^ = -Aiexp(-a;2) sinxi, B'^ = f3{x2) 


Q2 = Pi- ias 


Eq. (ISnD is 
shape invariant 


6. 


B^ = B^ = 0, 53 = /(x) 




decoupled 


7. 


B^ = B^ = 0, 53 = /(xi) 




decoupled 


8. 


B^ = B^ = 0, 53 = /(r) 




decoupled 



Here 



r='\Jxl + xl, 6' = arctan — , (28) 

hi-), /3(") and /(■) are arbitrary functions, and 5 = or 5 = 1. 

For B given in Items 1 and 4 to be single valued and Qi to generate finite rotations for 
arbitrary values of theta, parameter k must be an integer. 

The fields B presented in Items 2 and 3 of Table 1 are particular cases of the field presented 
in Item 1. We specify these particular cases since the corresponding vector B is divergent-free, 
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i.e., satisfies the condition 

V ■ B = (29) 

which is not vahd for more general fields given in Item 1. This condition is necessary if we 
interpret B as a vector of magnetic field strength. Analogously, the divergent-free external field 
presented in Item 5 is a particular case of the more general field given in Item 4. 

In the cases enumerated in Items 6-8 Hamiltonian is reduced to a direct sum of 
Hamiltonians with scalar potentials, and the corresponding equation ([T]) is reduced to the 
system of two decoupled equations. This fact is indicated in the fourth column of Table 1. 

Thus we find Hamiltonians ([3]) whose integrals of motion are generators of Lie symmetries. 
One more external field for which ([3]) admits Lie symmetries (but also a more general symmetry 
equivalent to with Ci = C2 = C4 = 0, C3 7^ 0) is presented in Item 1 of Table 2. 

Let us note that all obtained Hamiltonians admit a straightforward generalization which 
keeps its symmetries. Namely, in addition to (|3]), we can consider a more general Hamiltonian 

H = -V^ + a ■B + ujB'^ (30) 

where u is an additional coupling constant. 

Setting in ( l30l) w = we obviously come to the previous Hamiltonians ([3]). However, for 
all external fields presented in Table 1 the generalized Hamiltonians (J30l) admit the same Lie 
symmetries as Hamiltonian ([3]). 

One more generalization of Hamiltonian ([3]) which keeps its symmetries can be written as: 

H = -V"^ + (T -B + V (31) 

where V = V{r) and V = V{x2) for the cases presented in Items 1, 2, 8 and 3, 4, 7 correspond- 
ingly. 

V. NON-LIE INTEGRALS OF MOTION 

In this section we present integrals of motion which are not generators of Lie symmetries. 
Doing this we restrict ourselves to the first order integrals of motion whose general form is given 
by equation (j4]) where at least one of coefficients A"'', a = 1,2,3, is nonzero. Such integrals 
of motion can be represented in the forms given by equations ( l2Tll -( 1Ml) . Moreover, the latest 
equation should be supplemented by the condition + 7^ 0. 
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To evaluate the corresponding external fields B and functions fi*^ it is necessary to solve 
equations ( |271) where parameters a, b, ci, C2, C3, C4, di, d2 should satisfy one of the following sets 
of conditions: 

ab ^0, C3 = C4 = 0, (32) 
a = 0, by^O, di = d2 = 0, (33) 
a 7^ 0, 6 = 0, C3 = C4 = di = d2 = 0, (34) 
a = b = 0, cl + cl^O. (35) 

Moreover, conditions (152]) . ([2SD, (IM]) and (15^ correspond to symmetries (155]) . (12^ and 

(12^ respectively. In addition, equivalence transformations (ITSl) can be used to simplify the 
form of Q with fixed A'^'* and fi". 

Solving equations (!27|) under conditions (l32l) - (!35l) we obtain the list of external fields and 
the related integrals of motion which is presented in Table 2. 



Table 2. External fields and higher symmetries 



No 


External field 


Symmetry operators 


Comments 


1. 


= /icosxi, = /isinxi, 
B'^ = iy 


Q2 = Pi — lo's; -P25 

g3 = ^3(Pi-^^) 

— /i(cri COS Si + (72 sin Xi) 


supersymmetric , 
see (112]), 013]) 


2. 


j^l fiksm(k6) 

_ Mfccos(fce) ^3 _ ku 


Ql = L + |(T3, 
<54 = 0"3(<5l + '^) 

—fi{<Ji sm{k9) — (T2 cos(/i;6')) 


conformal, see (144]): 
single valued 
for integer/c 


3. 


fil - , 

52 - , 53 - 

2^i/2-^2,,2 2 


Ql = L + I0-3, 

Q5 = GiPi + 0-2^2 - f (criX2 

-o"2a;i) - |o-3a/^^ ^ /^^^^ 


supersymmetric , 
see 010]), dH]) 


4. 


r ' 

^2 ^ 53 ^ _^(^^y 


Qe = 0-1^1 + o-2^'2 + ^^{<yzQl 
+criX2(p - 0-2X1V2) + (7s{(p + v) 


supersymmetric , 
see (137]), (]38]), (I39]) 



Here = ip{r) is a solution of the following algebraic equation 

(/xV^ + 1) ^ = c, (36) 
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fi, u, k and c are real parameters. In particular, for z/ = 0, c = and c = 0, u = —Au we 
obtain the following versions of the field presented in Item 4: 

, u sin 9 sinh p , w cos ^ sinh p „3 

is = o , ij" = T , -D = Q — 

cosh p cosh p cosh p 

and 

w sin 9 sinh p ^2 ^ cos 6* sinh p ^3 a; 

cosh'^p ' cosh'^p ' cosh^p 

correspondingly, where we denote r = sinh p/p. 

Thus we have found the complete set of integrable planar models of neutral particles with 
Pauli interaction. The Hamiltonians of these models are given by equation (jS]) where B is the 
vector of external field whose components are presented in Tables 1 and 2. 



VI. ALGEBRAS AND SUPERALGEBRAS OF SYMMETRY OPERATORS 

Symmetry operators presented in any item of Table 1 and Table 2 commute each other. In 
other words, the presented sets of integrals of motion together together with the corresponding 
Hamiltonians form bases of Abelian Lie algebras. 

In addition, integrals of motion collected in Table 2 form interesting superalgebraic struc- 
tures. Namely, operators presented in Item 4 satisfy the following relations: 

Ql='H, [Q6,H] = 0, (37) 

[QuH] = [Qi,Q2] = 0. (38) 

where 

n = H +{pQi + uf + c (39) 

and H is the corresponding Hamiltonian 1^. In other words, operators Qi, Q2 and H form 
a basis of the Lie superalgebra whose odd and even basis elements are Q2 and < QijTi > 
correspondingly. Relations f l37|) specify = 1 SUSY. 

Relations (1371) and (138!) can be effectively used to find eigenvectors and eigenvalues of the 
corresponding Hamiltonian ([3]). Indeed, the commuting hermitian operators Qi, Q2 and H 
have common eigenvectors. To find eigenvectors for the first order differential operators Qi 
and Q2 is much more easier than for the Hamiltonian which is a differential operator of the 
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second order. In addition, relation f p7|) makes it possible to find eigenvalues of H algebraically 
starting with eigenvalues for Qi and Q2- 

Operators presented in Item 3 of Table 2 together with the corresponding Hamiltonian (|3]) 
also form a basis of superalgebra since the following relations are satisfied: 

(40) 

[Qi,H] = [Qi,Q5]=0 

where 

'H = H + fiQi. (41) 

In addition to their mutual commutativity the symmetry operators presented in Item 1 of 
Table 2 satisfy the following quadratic relations 

Ql = H + 2uQ2 + iy'', (42) 

Q3-2J =(Q2 + ^)'. (43) 

Thus the corresponding models also admit = 1 SUSY and can be effectively integrated. 
Integrals of motion presented in Item 2 of Table 2 satisfy the following algebraic relations: 

Ql = Ql + 2uQ, +fi^ + iy^ (44) 

which can be used to find eigenvalues of Q2 using (well known) eigenvalues of Qi. Moreover, 
these integrals of motion together with operators D = XiPi + X2P2, K = r^/2 and the corre- 
sponding Hamiltonian H form a basis of the five- dimensional Lie algebra since the following 
commutation relations are satisfied: 

[H, D] = -2iH, [K, D] = 2iK, [K, H] = iD (45) 

while all the other commutators are trivial. 

Relations (H5ll characterize conformal algebra so(l,2), thus we deal with a model of conformal 
quantum mechanics ( for definitions see, e.g., 22|). 



VII. EXACT SOLUTIONS 



In this section we use symmetries of the models found above to construct their exact so- 
lutions. Thanks to the presence of arbitrary parameters the number of qualitatively different 
models is too large to be considered in one paper, and it is the reason why we restrict ourselves 
only to three particular examples. 
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A. Neutron in periodic magnetic field 



Let us start with the relatively simple Hamiltonian whith the components of magnetic 
field B given in Item 1 of Table 2. Since this magnetic field depends on one spatial variable, it 
is reasonable to restrict ourselves to the one- dimensional eigenvalue problem 

Hip = ( + /^(^i cosy + (T2 siny) + uasj ip = Eip (46) 

where we denote xi = y. 

Hamiltonian in fH6l) includes periodic potential, thus this equation is a certain analog of the 
Bloch problem for electron. However, the considered eigenvalue problem is related to a neutral 
particle and includes a spin dependent potential. 

Hamiltonian (146 p commutes with operators Q2 and Q3. Thus we can search for common 
eigenf unctions for Q2, Q3 and H. 

The eigenvalue problem for the first order differential operator Q3: 

Qsi^k = (0-3 (A - z/) - fi{ai cos(y) + (72 sm{y))) ip = kipk (47) 
is easily solvable. The general solution for equation fl47p is the two component function 




= r I (4j 

where 



(l)i{y) = exp ( '^'^''^'^^y ) ((Cik_ + C2Xk) cos(Afcy) + (CiA^ + Cafc.) sm{Xky)), 



'ii2u l)y\ 

^2(2/) = -/iexp ( ) (Cicos(Afc?/) + C2sin(Afcy)). 



Here k_ = k — ^, A^ = a/A?~^^7'^, and the latest quantity must be real if we ask for solutions 
whose norm does not turn to infinity with growing y. In other words, admissible values of k 
are restricted by the condition 

V 



k-l) > fi^. (50) 



Let us present the admissible values of k more explicitly. First we note that up to the unitary 
transformation H — > a-^Ha^i we can restrict ourselves to /i > 0. Then, in accordance with fISU]) 
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there are two possibilities: 



k>n + l or k<-ix + ^. (51) 



Since operator Q2 commutes with H, eigenfunctions (HHj) . ( H9|) solve also equation ( H6|) . The 
corresponding eigenvalues E are easily calculated using algebraic relation = H: 

E = k\ (52) 

In accordance with fl5T]) . fl52l) the admissible values of E are restricted by the following 
conditions: 

E>\fi--j if ^>2' ^<2~^' 

E>0 if 0</i<^, k<^-fi, (53) 

The probability density corresponding to solutions f l49|) . i.e., 

0101 + 020; = (C^? + C|)P + 2A;_A,CiC2 + ^(^2 - C|)^2 cos(2A,i/) (54) 



is a periodic function. However its period can differ from the shift which keeps equation fl46|) 
invariant, i.e., from 2tt. Such situation looks rather nonphysically, since the probability density 
calculated in a fixed frame of reference can differ from the density calculated in the equivalent 
frame of references shifted by 27r, in spite of that the equation in these frames has exactly 
the same form. 

There are two ways to obtain solutions whose amplitude is a periodic function with the 
period 27r. First it is possible to choose in (149|1 



Ci = ±C2 = ^ = (55) 

2v/vrA;_(A;_ ± Afc) 



and obtain solutions normalized at any invariance interval [y, y + 27r]: 



In this case there are no restrictions on eigenvalues E additional to ( 1531] . 

14 



The second way is to impose the following condition on the spectral parameter k: 

k = ^(5V«2 + V + 1) (56) 

where n = 0,l,2,..., e = ±l. In this case the energy levels fl52p are discrete and there are 
two branches: and where 

E± = hji^ + V ± 2Vn2 + V + !)• (57) 

The corresponding eigenfunctions can be obtained from (HSj) . (H^ changing — )■ | and using 
expression ( 156|) for fc. To obtain normalized solutions, arbitrary constants C\ and C2 should be 
restricted by the following condition: 

2 



TT 



Thus like in the Bloch problem for electron 23| the energies of neutron moving in the periodic 
field represented in Item 1 of Table 2 can be continuous and have a band structure, see fl53l) . 
In addition, there are solutions with discrete spectrum (!57j) . 



B. Rotationally invariant system 

The next system which we consider includes the following Hamiltonian: 

= -V^ + ^{axX2 - o^xx) + - (58) 

where /i and a are real parameters. 

If a = Hamiltonian (135]) coincides with operator ([3D where and B"^ are components of 
the external field presented in Item 2 of Table 2 with k = \ and i/ = 0. For nonzero a this is a 
Hamiltonian of type fl3Tl) . 

The additional term - which we include to obtain a more general model does not break the 
commutativity of the Hamiltonian with operators Q\ and Qt^. Using this fact and taking into 
account relations (jSj) we can expand solutions of the eigenvalue problem ([T]) for Hamiltonian 
( l58l) via eigenvectors of these operators satisfying 



( CT3L - -(criX2 - cr2Xi) + - j ^fc,^ = + /i^^fc,^, e = ±1. 
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Solutions of equations fl5^ can be represented in the following form 

_ Cfe.0(r) f expjtjk - 1)6) {k + gy/P + f,^y ^^^^ 
iexp(z(/c + i)6')/i 



where Cke are integration constants and polar coordinates (125]) are used. Substituting 
( l60l) and ( l28l) into ([T]) we obtain the following ordinary differential equation for radial function 
(f) = 0(r): 



2 



Let a = then equation (I6T1) with the boundary condition 

= if r = (62) 

defines the eigenvalue problem of one dimensional conformal quantum mechanics. Its solutions 
can be expressed as follows: 

E = / > 0, (j) = ^J^{pr) (63) 
where Jy{pr) is the Bessel function with 



V = ^^1 + 4P-4£VP + /x^. (64) 
If e = 1 then admissible values of parameters fi and k are constrained by the following relation: 

fc^-^j >/^^ if e = l. (65) 

Let a > then solutions of equation fl^Tl) can be found using its analogy with the radial 
equation for Hydrogen atom. The eigenvalues E which correspond to square integrable solutions 
vanishing at the singularity point r = look as follows: 

4 (n + I + z/) 

where u is parameter given in ( l64l) . n is a natural number and k satisfies conditions ( l59l) and 
0651) . The corresponding eigenvector can be expressed via a linear combination of Whittaker 
functions M{a,b,x) and W{a,b,x): 

= CiM (^a, 6, ^) + C2W (a, 6, ^) (67) 
16 



where 

3 , 1 
a = n + u + -, b = u + -. 

4 4 

Notice that in contrast with the Hydrogen atom there is not a degeneration w.r.t. orbital 
quantum number. However, any energy level is infinitely degenerated since the correspond- 
ing eigenvector (|671) includes two integration constants, an only one of them can be fixed by 
normalizing the wave function. 



C. Shape invariant system 

Let us consider the eigenvalue problem for Hamiltonian (130|) where B is the magnetic field 
whose components are given in Item 4 of Table 1 where /s = 0: 

Hip = (-V^ + A(l - 2K)exp(-X2)(o-icosa;i - sin xi) + A^exp(-2x2)) ip = Eip. (68) 

Here A is the integrated coupling constant and new parameter k is introduces such that /i = 
A(l - 2k). 

Hamiltonian H in (168!) admits integral of motion Q2 = Pi — Thus it is possible to expand 
solutions of f l68|) via eigenvectors of Q2 which look as follows: 

y exp{i{p - |)xi)^(x2) J 

and satisfy the condition Qipp = pipp. 

Substituting f l69|) into f l68|) we come to the following equation: 



+ K $ = (70) 



where we denote 



2 1 ^ 



y = X2, e = E-p'--, $=|';| (71) 
and 

K = A^ exp(-2y) - X{2k - 1) exp{-y)ai - pa^. (72) 
We can restrict ourselves to non-positive p, then solutions for p > could be obtained changing 

$ (Ti$. 
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Let us consider the eigenvalue problem (I7U1) with the following conditions: 

POO 

y>0, $(0)=0, / $*(i/)$(|/)(i|/<cx). 

^0 



Potential (172]) belongs to the list of shape invariant matrix potentials found in 20|], see 
equation (5.12) therein. It can be represented in the form V = — W'^ + where 

2 

W^, = -K + X exp{-y)ai - ^^"3, = + (73) 

The shape invariance means ^ that the superpartner potential V^"*" = + + is equal 
to the initial potential with shifted parameter k up to a constant term. And it is the case for 
potential (1721) since 



Using the shape invariance it is possible to integrate equation (|68|) in a simple and straight- 
forward way with using tools of SUSY quantum mechanics [20!]. The eigenvalues e and the 



corresponding state vectors are enumerated by natural numbers n = 0,1, . . . The ground state 



vector $o(K,y) = I I should solve the equation 



a-%{K, y)=(J- + W^^ y) = 0, (74) 

thus 

(^0 = z'^-'^K.+^iz), ^0 = z^-^K,{z) (75) 

where Ky{z) is modified Bessel function, = ^ — | and z = Aexp(— y). 

Solutions which correspond to n^^ exited state can be calculated using the following relation: 

$„(k, y) = a+at+i ■ ■ ■ a^+n-M^^ + v) (76) 
where a'^ = — ^ + W^iy)- Finally, the corresponding values of spectral parameter e have the 



following form 20| 



e = -N'-^^E = p^-N^-^ + ] (77) 



4iV2 ^ Am A 

where N = k + n and n is a natural number. 
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In accordance with fl77j) eigenvalues E are invariant w.r.t. the change p — > —p, thus eigen- 
vectors corresponding to a chosen n are hnear combinations of functions (l69l) : 



^ /exp(i(p+ i)xi)¥?„(x2)\ f exp(i(-p + i)xi)^„(x2) 

yexp(z(p - |)xi)^„(x2) y yexp(i(-p - |)a;i)v3„(x2) 

Let C1C2 7^ then, in order the norm of this function be invariant w.r.t. the shifts Xi — )■ 
Xi + 27r hke the Hamiltonian f l68|l . it is necessary to impose the following condition: 

2m + 1 
p= , m = 0, 1,... 

For solutions ( 1781) with Ci = or C2 = the spectral parameter p can be quantized by imposing 
the periodic boundary condition with an arbitrary period. 



VIII. DISCUSSION 



In the present paper planar Schrodinger-Pauli equations for neutral particles, which admit 
first order constants of motion, are classified. The collection of such equations appears to be 
rather rich and interesting. In particular, it includes supersymmetric systems belonging to 

= 1 SUSY quantum mechanics, the system with shape invariant Hamiltonian (!68|) and the 
system with Hamiltonian f l58|) (were a = 0), which is conformally invariant. 

Any Hamiltonian ([3]) with the external fields presented in Table 2 admits two integrals of 

71 n 

motion, and following pj| -p^ we can call the related 2d systems superintegrable. Notice that 
this terminology is rather conventional, since the additional (spin) degree of freedom is ignored. 
To fix this degree of freedom we need an additional involutive integral of motion like matrix a^, 
which can extend the possible number of symmetries. Maybe it is more natural to say that a 2d 
system with spin 1/2 is "superintegrable" (integrable) if it has at least three (two) independent 
integrals of motion. These speculations can be justified by the example given in Item 6 of 
Table 1. The corresponding 2d system admits a constant of motion, but it is not convenient 
to call it "integrable" since the related Hamiltonian is a direct sum of two 2d " nonintegrable" 
Hamiltonians which do not admit constants of motion provided is an arbitrary function of 
Xi and X2- 

Quantum mechanical systems with a sufficiently large number of constants of motion are 
usually exactly solvable, and it is the case for the models classified in the above. We restrict 
ourselves to solving three of them in section 7. To solve the model whose Hamiltonian is given 

19 



by equation (1^ we use its shape invariance, and this is the second direct apphcation of results 
of paper {20 1 where matrix superpotentials were classified, to a d-dimension models with d> 1. 
The first application of these results to planar systems with arbitrary spin can be found in 
paper 

A physically interesting subclass of the classified systems includes Hamiltonians ( H6|) and 
(jnSl) whose effective potentials are periodic functions. Such potentials simulate interaction of 
neutron with a crystal lattice. 

Fundamental results concerning the motion of electron in periodic electromagnetic field 
was formulated long time ago by F. Bloch [23]. Our analysis of solutions f H9|) shows that the 
neutron interacting with a periodic magnetic field can have both continuous and discrete energy 
spectrum. 

The presented classification of planar Hamiltonians admitting first order integrals of 
motion is complete. Nevertheless it can be considered as an intermediate outcome. First it is 
interesting to study integrable planar systems with higher order integrals of motion. An example 
of such system is the Pronko-Stroganov model [l^ which admits symmetry operators of second 
order. Then, our analysis can be extended to systems of type fl30|) with arbitrary B and V. 
Finally, the 3d sup erintegr able systems for neutrons are also waiting for their classification. In 



other words, superintegrable systems with spin whose investigation was started with paper 
belong to a promising research field. 
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